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ABSTRACT 

We study a non-dissipative hydrodynamical mechanism that can stabilize the spin of the ac- 
cretor in an ultra-compact double white dwarf binary. This novel synchronization mechanism 
relies on a nonlinear coupling between tides and the uniform (or rigid) rotation mode, which 
spins down the background star. The essential physics of the synchronization mechanism is 
summarized as follows. As the compact binary coalesces due to gravitational wave emission, 
the largest star eventually fills its Roche lobe and accretion starts. The accretor then spins 
up due to infalling material and eventually reaches a spin frequency where a normal mode 
of the star is resonantly driven by the gravitational tidal field of the companion. If the res- 
onating mode satisfies a set of specific criteria, which we elucidate in this paper, it exchanges 
angular momentum with the background star at a rate such that the spin of the accretor locks 
at this resonant frequency, even though accretion is ongoing. Some of the accreted angular 
momentum that would otherwise spin up the accretor is fed back to the orbit through this 
resonant tidal interaction. In this paper we solve analytically a simple dynamical system that 
captures the essential features of this mechanism. Our analytical study allows us to identify 
two candidate Rossby modes that may stabilize the spin of an accreting white dwarf in an 
ultra-compact binary. These two modes are the / = 4, m = 2 and Z = 5, m = 3 CFS unsta- 
ble hybrid r-modes, which, for an incompressible equation of state, stabilize the spin of the 
accretor at frequency 2.6a;orb and 1.54 Worb respectively, where o^orb is the binary's orbital 
frequency. For an n = 3/2 polytrope, the accretor's spin frequency is stabilized at 2.13 Woib 
and 1.41 cjoib respectively. Since the stabilization mechanism relies on continuously driving 
a mode at resonance, its lifetime is limited since eventually the mode amplitude saturates due 
to non-linear mode-mode coupling. Rough estimates of the lifetime of the effect lie from a 
few orbits to possibly millions of years. We argue that one must include this hydrodynamical 
stabilization effect to understand stability and survival rate of ultra-compact binaries, which 
is relevant in predicting the galactic white dwarf gravitational background that LISA will 
observe. 

Key words: white dwarfs - binaries: close - novae, cataclysmic variables - gravitational 
waves - stellar dynamics. 



1 MOTIVATION, MAIN RESULTS AND OUTLINE 

The collection of ~ 100 to 200 million double white dwarf bi- 
naries (WD-WD) populating the Galaxy generates an important 
gravitational wave background that the planned LISA mission will 
be sensitive to. Predicting the properties of this gravitational wave 
background requires accurate population synthesis models, which 
in turn necessitate precise understanding of dynamics of WD-WD 
binaries, which we now briefly summarize. 

Upon formation a WD-WD binary will have its constituents 
separated widely enough so that no mass transfer is occurring. Dur- 
ing this detached phase the orbital motion of the binary is well ap- 
proximated by that of point-particles in Newtonian gravity, sup- 
plemented by leading-order tidal coupling and dissipative contri- 



butions due to gravitational wave emission, computed from the 
Burke-Thorne potential. As the binary coalesces, the larger white 
dwarf will eventually fill its Roche lobe, signaling the onset of mass 
transfer. A sizable fraction of galactic WD-WD binaries are com- 
pact enough so that this mass transfer phase will be reached within 
a Hubble time. During mass transfer the binary's dynamics is sub- 
stantially more complicated than during the detached phase, since 
the accretion rate and orbital parameters are coupled to one another. 
In particular the question of whether or not WD-WD binaries sur- 
vive mass transfer for a long time, as opposed to merging shortly af- 
ter mass transfer begins, is not fully understood at the present time. 
If indeed such compact binaries are stable, then they will be among 
the strongest persistent sources contributing to the galactic gravita- 
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tional wave background for LISA. There are currently at least two 
observed ultra- compact WD- WD candidate system s, namely RX 
J0806.3+1527 ("StrohmaveJ dlOOSh ) and V407 Vul iRamsav et al\ 
( l2005h : IStrohm aver (2004i)). 

The problem of stability of accreting white dwarf binaries ha s 
been investigated in detail bv lMarsh. Nelemans & Steeghj j2004h . 
They show that for mass ratios lying between guaranteed stabil- 
ity and guaranteed instability, the stability of accreting WD-WD 
binaries is closely related to the strength of the dissipative synchro- 
nization torque Xiuss that couples the accretor's spin and the orbit, 
which is assumed to be of the form 



(3) 



(1) 



where rs is the synchronization timescale, / the accretor's moment 
of inertia, Q the spin frequency of the accretor and uiovh the bi- 
nary's orbital frequency. Their result is that for accreting WD-WD 
binaries to be stable the synchronization timescale must be < 10^ 
years. 

The essential contribution of the present paper is to highlight a 
non-dissipative hydrodynamical mechanism that may affect signifi- 
cantly the evolution of the spin freque ncy of the accretor, and there- 
fore im pact the stability analysis of 'M arsh. Nelemans & Steeghj 
Specifically we show that resonant tidal excitation of gen- 
eralized r-modesQ (or Rossby modes) in the accretor may stabi- 
lize its spin at a given resonant frequency of order orbital fre- 
quency. During this typ e of hydrodynamical resonance locking 
jWitte & Savonii^ j 19991) ') phase, some of the angular momentum 
carried by accreted material that would contribute to spin up the 
white dwarf is instead fed back into the orbit, which contributes to 
make the binary more stable to mass transfer. 



1.1 The tidal stabilization mechanism 

The physical foundation of the stabilization mechanism we ana- 
lyze is the coupling of the uniform rotation mode with other stellar 
modes excited by the companion's tidal field. Formally this cou- 
pling c omes out natu r ally fr om second order stellar perturbation 
theory jSchenk et all ( l2002h ). but may be understood schemati- 
cally as follows. Consider a perturbation of a uniformly rotating 
star characterized by a fluid Lagrangian displacement vector field 
^. The total angular momentum of the perturbed star may be com- 
puted as an expansion in ^, which can be formally written as 



Jstar = Jo + Mi] + J2[i,i], 



(2) 



where Jo is the total angular momentum of the unperturbed star. 
Perturbation theory of uniformly rotatin g stars yields that J i de- 
pends only on the uniform rotation mode jSchenk et ai\ ( l2002h ). all 
other stellar modes contributing only in J2. This implies that the 
first two terms in the right-hand side of l|2j can be meaningfully 
combined into a term of the form Jfi, where / is the moment of 
inertia of a star rotating uniformly at spin frequency Q, since the 
angular momentum Jo + Ji is carried entirely by the uniform rota- 
tion modfl Taking a time derivative of (|2]l then yields 



^ Generalized r-modes are predominantly toroidal perturbations whose 
restoring force is provided by the Coriolis force, which makes them dy- 
namically unimportant in non-rotating stars. 

^ However in situations where the unpeiturbed star has significant differ- 
ential rotation, it is not clear if this result canies over as the formalism of 
ISchenk et al.\ j2002l) has not been yet extended to unperturbed stars with 
differential rotation. 



In binary systems one may compute the contribution to the exter- 
nal torque on the star due to non-dissipative tidal coupling as an 
expansion in ^ as well, using either the binary's equation of motion 
as done in Appendix [D] or equivalently the tidal interaction poten- 
tial. The total torque on the accretor is then the sum of this tidal 
torque — Jtidai and the accretion-induced torque — Jacjfl leading 
to the following angular momentum conservation equation 



^^{m) + jaK,^] = -(Jtidai + Jacc). 



(4) 



Expanding the lagrangian displacement i{x,t) in terms of normal 
modes {jj,(a;jf|as 



i{x,t) = 2^Ca{t)i„{x) + CO., 



(5) 



one finds that — Jtidai does not equal Jaf^,^], the difference be- 
tween both terms describing how the uniform rotation mode cou- 
ples to other stellar modes. This coupling however should not de- 
pend on whether the perturbed star is part of a binary system or 
not. Thus there should exist a derivation of this coupling contained 
entirely within the framewor k of stellar pe r turbat ion theory. This 
derivation is indeed given in ISchenkefa/.! i l2002h and yields the 
same answer as Eq.([4l( when specialized to a binary system. 

Now in studies of accreting binary systems it is generally as- 
sumed that the difference between the terms — Jtidai and J2 in ^ is 
negligible compared to the accretion-induced torque Jacc . In situa- 
tions where the tidal response of the accretor's normal modes to the 
gravitational field of its companion lies away from any resonance, 
this yields a very good approximation. However as the accretor 
accumulates material and spins up, it becomes possible to sweep 
through tidal resonances, in which case the back-reaction of the 
resonating mode on spin frequency evolution may not be neglected. 
The mode amplitude grows on a short timescale and the difference 
between J2 and — Jtidai may become comparable in magnitude to 
Jacc . As we show in this paper, for a simple model where a single 
mode dominates the tidal response of the star, Eq.(|4ll assumes the 
form 



at I 



(6) 



where Ca is the amplitude of the resonating mode and I'a is a pa- 
rameter describing the coupling between the resonating mode and 
uniform rotation. This parameter typically scales as i>ci ~ <^a, 
being the mode eigenfrequency. In the context of a binary system 
one may also interpret the second term in the right-hand side of ((SJ 
as the difference between the external tidal torque due to the com- 
panion and the angular momentum carried by the resonant mode. If 



Here by "non-dissipative tidal coupling" we simply mean the coupling 
between the mass multipole moments of the accreting star induced by the 
fluid perturbation and the Newtonian gravitational tidal field of the com- 
panion. In this paper we omit dissipative tidal coupling for the puipose of 
computing the tidal response of the accretor We do however mention it 
when discussing stability issues later in the introduction. 
* We introduce minus signs here to follow the convention used in the body 
of the paper, which defines Jtidai and Jacc as rates of change of orbital 
angular momentum due to tidal interactions and accretion respectively. This 
implies in particulai' that Jacc must be negative. 

^ In this work the mode eigenfunctions have dimensions of length, so the 
mode amplitude coefficients Ca ai'e dimensionless. The normalization con- 
dition on the eigenfunctions we use is detailed in Appendix Icl 
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Va is positive and large enough, it may then be possible that when 
mode amplitude Ca grows rapidly due to resonant excitation, the 
accretion-induced torque Jacc is balanced by the hydrodynamical 
back-reaction term, resulting in a spin that is nearly constant in 
time. The value of this quasi-static spin frequency is determined by 
the eigenfrequency of the resonant mode. 

Let us again emphasize that this spin stabilization mechanism 
is entirely hydrodynamical since it relies solely on tidal excitation 
of normal modes of a perfect fluid. No internal damping is needed 
to stabilize the spin in this scenario. However, as opposed to dis- 
sipative synchronization, the lifetime of this mechanism is limited 
since it requires continuously driving a mode near resonance. Even- 
tually the mode amplitude will grow large enough so that non-linear 
mode-mode coupling cannot be neglected. Generically mode-mode 
coupling sets a saturation amplitude beyond which it is not possi- 
ble to drive the mode effectively. The lifetime of our stabilization 
mechanism therefore depends crucially on this saturation amplitude 
which, for resonant tidal excitation of modes, is not well known. A 
computation of this saturation amplitude, as well as understanding 
the evolution of the binary when the mode saturates, likely requires 
simulating a large network of coupled modes following the work 
of ISrink, Teukolskv & Wasserman (2004, 2005), who studied the 
problem of growth of r-modes unstable to gravitational radiation 
reaction. As this is a very complicated dynamical system, we shall 
remain cautious and not speculate any further on the fate of the 
binary after saturation. 



1.2 Overview of the results 

The central result of this paper is the proof that there exists modes 
in a rotating star whose resonant tidal excitation can stabilize the 
spin of the accretor in ultra-compact binary white dwarf systems. 
This result is derived following a number of steps. First we solve 
analytically the coupled system of equations describing the time 
evolution of the spin frequency of the accretor Q, and the mode 
amplitude Ca ■ This system may be written as 



- d 2 



- J- a exp 



Q.{t')dt' 



(7) 



(8) 



where J^a is the overlap of the mode eigenfunction with the ex- 
ternal tidal field and (f>orh is the binary's orbital phase. To solve 
system l|7j-l[8j we make the following three simplifying approxi- 
mations: (i) we assume the accretion-induced torque is constant, 
(ii) we neglect the back-reaction of the modes on orbital evolution 
for the purpose of solving l[8} and (iii) we assume the accretor is 
slowly rotating, i.e. f7 <C M/R^. [Throughout we employ geo- 
metric units G = c = 1.] In the case of a mode which stabilizes the 
spin frequency, the analytic solution for the spin frequency during 
resonance locking may be written as 



Aq| Jaccj {t — to) 



(9) 



where fires is the spin frequency at which mode a resonates, is 
a dimensionless number of order unity and time to is defined as the 
time when the mode enters resonance regime. 

Once system ([7} is solved, we investigate whether or not solu- 
tion (|9} is physically realized by some set of normal modes of the 
accretor. We show that solution ^ is only valid for resonantly ex- 



cited modes whose back-reaction parameter and tidal coupling 
strength J^a satisfy the following bound 



J a< 



3/2 



(10) 



In order to stabilize the spin of the accretor at frequency of order 
orbital frequency, one needs to focus attention to modes with low 
corotating-frame eigenfrequencies, as the stabilization frequency 
Sires is related to orbital frequency cjorb by 



S^rcs — 



(11) 



The ideal candid ate modes ar e generalized (or hybrid) r-modes 
(Rossby modes) tervani ( Il889h : iLindblom & ipserl ( Il999l) ), since 
their eigenfrequencies scale as uia ~ f2. Therefore they are very 
low frequency modes in slowly rotating stars, compared to, say, / 
modes. We then show that in ultra-compact binary white dwarf sys- 
tems the following generalized r-modes are potential candidates for 
stabilizing the spin of the accreting star: 

(i) The I = 4,m — 2 r-mode with eigenfrequency uji2 = 
— 1.232 This mode stabilizes the spin at frequency fires, 42 = 

2.60 Uorb- 

(ii) The I — 5, m — 3 r-mode with eigenfrequency LJ53 — 
— 1.053 fi. This mode stabilizes the spin at frequency fires, 53 ~ 

1.54 Uorb. 

There exist other candidate modes with polar quantum number / = 
6 but they satisfy bound l llOt only marginally. 

As mentioned previously Eq.l|9]l is valid only up to a cer- 
tain time tmax when the resonant mode amplitude becomes large 
enough so that non-linear hydrodynamical effects (mode-mode 
coupling) saturates its growth. In that regime the simple model we 
analyze here does not describe correctly the tidal response of the 
star. As mentioned previously the fate of binary in the saturation 
regime remains an open problem. Denoting the saturation mode 
amplitude as Csat , the lifetim^fl Tm of the hydrodynamical stabiliz- 
ing mechanism is then found to be 



Tiif 



Csat 



J acc 



(12) 



The current uncertainty in the saturation amplitude and different 
choices of orbital parameters together give estimates for the life- 
time of the mechanism ranging from a few orbits to hundreds, pos- 
sibly millions of years. We also investigate how our results change 
when the stellar model is an n = 3/2 polytrope. In that case we 
find that backreaction of the ? = 5, m = 3 CPS unstable r-mode 
may or may not be strong enough to stabilize the spin, depending 
on accretion rate and compactness of the binary. If it is not strong 
enough, then the / = 4, m = 2 CPS unstable mode will be able to 
stabilize the spin in the case of an n = 3/2 polytrope. 

1.3 Implications of the stabilization mechanism on binary 
stability and observations 

Before going ahead with the detailed derivation of the re- 
sults presented above, we will discuss briefly some observa- 
tional consequences of our tidal stabilization mechanism. More 

By lifetime we here mean the time interval during which the tidal re- 
sponse of star and the distribution of angular momentum is correctly de- 
scribed by the model studied in this paper. 
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specifically we take a look at h ow the stability analysis of 
iMarsh. Nelemans & Steeghsl ( l2004h is modified when the binary is 
in the resonance locking regime. During resonance locking the spin 
frequency of the accretor is essentially constant, so Eq.l|6) yields 



I I 2 acc 

dt 



(13) 



The evolution equation of the orbital angular momentum is given 
by 



^orb — JgW + Jacc + Jtidal.diss ' 



, d I |2 



(14) 



where Jqw is the angular momentum loss due to gravitational 
waves, Jtidai.diss is the dissipative tidal torque usually modeled 
as in Eq.([T) and the last term of the right-hand side is due to 
modal tidal coupling [see AppendixlD). The number be, sometimes 
called wave action, is mode energy at unit amplitude divided by 
corotating-frame mode frequency. Substituting l ll3t into il4\ yields 



Jorh = JgW + (1 — Xa) Jacc + Jtidai.diss, 



where 



(15) 



(16) 



For the two r-modes mentioned previously, ha is negative and 
therefore Xa is positive, and is interpreted as the fraction of ac- 
creted angular momentum that is fed back to the orbit by the reso- 
nant tidal interaction. It turns out that for both modes Xa is approx- 
imately equal to 0.5. The binary is more stable to mass transfer 
during the resonance locking regime, as the coefficient in front of a 
term driving an instability is reduced. Stability of the binary under 
mass transfer is governed by the evolution of the Roche lobe over- 
fill parameter A = R2~ Rl, where R2 is the donor radius and Rl 
is the donor Roche lobe size. The equation governing the evolution 
of the over-fill parameter during resonance locking is given by 



1 dA 
2R^~dt 



(C2-C.L) 



M2 
Mo 



= J ^ (fi - '^orb) +cr—, (17) 

where q is the mass ratio, A/2 is the donor mass and a is defined as 

(18) 



= 1 + i(C2 - Cr J - g - (1 - Xa)y/{l+q)Rh/ 



The quantity Rh is the radius of the orbit around the accretor that 
has the same specific angular momentum as the accreted mass, 
and can be conveniently ap proximated by the following expression 
jVerbunt & RappaportI l fT988l) ) 

— = 0.088 -0.049 log <3 + 0.115 log^ (7 + 0.020 log^g, (19) 
a 

which is valid for the range 0.001 < q < 1. For mass ratios ranging 
from 0.1 to 1, we have Rh/a ~ 0.15 in order of magnitude. Lastly 
we have 



C2 



d log R2 
d log M2 ' 

d\og{RL/a) 
d log M2 ' 



(20) 



(21) 



with typical values —0.6 < C2 < —0.3 and — 1/3. The binary 
will be unstable if mass transfer AI2 < increases the over-fill. 



Figure 1. This figure shows the boundary of the regimes of guaranteed sta- 
bility in two cases: (solid line) when no modes are driven, which is the 
criterion of Marsh, Nelemans and Steeghs (2004), and (dashed line) during 
resonance locking, which is Ea. j22t with Xa = 0.5. The binary is stable 
to mass transfer in the region of the graph below the line corresponding 
to the appropriate regime. Clearly resonance locking increases significantly 
the parameter space over which the binary can undergo stable equiUbrium 
mass transfer 



since a larger overfill implies in turn a larger mass transfer and thus 
the process runs away, leading to the destruction of the binary. 

Let us now look for quasi-equilibrium solutions to J17t . i.e. 
solutions to A = 0. The first term on the right-hand side of l ll7t 
is positive and the second one is negative, since during resonance 
locking the (constant) spin frequency is larger than the orbital fre- 
quency. However in the limit of very weak dissipative tidal cou- 
pling (ts ~ lO'^^ yr), which we assume in this paper, the second 
term can be neglected and the criterion for existence a stable equi- 
librium solution is simply ct > 0, which translates to 

q < 1 + ^{C2 - (rr.) - (1 - Xa)y/il + q)Rh/a. (22) 

Using Eggleton's zero-temperature m ass-radius relation and his ap- 
proximation for as both quoted by Mars h. Nelemans & Steeghsl 
(20.Q4|) ; we can compare stability criterion ll22t to the stability crite- 
rion of lMarsh. Nelemans & Steeghsl ( |2004|) [Eq.(3 1 ) of that paper] ; 
the result is shown in Figure[T] 

Now if a < 0, we have A > (still neglecting dissipative 
tidal coupling) and no stable equilibrium solution exists during res- 
onance locking. In that case, the accretion rate will keep increasing 
until it becomes too strong for hydrodynamical spin stabilization to 
operate, and the accretor can potentially be spun up through the res- 
onance without saturation of the mode amplitude. A precise mod- 
eling of the evolution of accretion rate during resonance locking is 
required to accurately quantify this effect, which we do not address 
in the present paper. 

It is also interesting to ask what the observational signature of 
the resonance locking regime would be, for example in the evolu- 
tion of orbital frequency. Before resonance occurs the effect of tidal 
coupling due to modes on orbital angular momentum evolution is 
negligible. The equation governing the rate of change of orbital an- 
gular momentum in this regime is thus Eq.lllSt, with Xc, taken to be 
zero. As resonance is approached, we show in section|4]and figure 
|2]that the accretor' s spin gets stabilized over a timescale of about 
a year. Once the spin is stabilized the orbital angular momentum 
evolves according to l IlSI l with Xa — 0.5 for the modes considered 
here. This implies that one should observe a sudden decrease SuJorh 
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in the time derivative of tlie orbital frequency ujorh at the onset of 
resonance locking regime, given by 



In this section we also identify relevant candidate modes for stabi- 
lization and perform the estimate of the mechanism's lifetime. 



5u)o 



oXa J tVacc 
"orb 
'i'orb 



<-'orb 



(23) 



which is negative since J^cc < 0. Since ilioTh is negative because of 
mass transfefl onset of resonance locking produces an increase in 
the magnitude of cjorb ■ Alternatively one may rewrite l l23t in terms 
of orbital period derivative Porh as 



SPoih — 



" orb 



6 X 10" 



IM2I 



lO-5M0yr-i 



Ah 

lOminy \^iW0 

-1 

miri yr , 



(24) 



where we have assumed Xa ~ 0.5 and the order of magnitude 
Rh ~ 0.15a to obtain the second line. Once the resonant mode 
saturates, the stabilization mechanism turns off and we may essen- 
tially reset Xa to zero in l ll5t . leading to a rapid increase in ti^orb, 
given by the negative of ( I23t . Unfortunately this observational sig- 
nature does not correspond to what is seen in RX J0806.3+1527 
and V407 Vul, where the orbital frequency is observed to be in- 
creasing with time at a rate consistent with orbital dynamics being 
entirely dominated by gravitational wave radiation reaction. Reso- 
nance locking due to tidal excitation of r-modes tends to push the 
stars apart even further than mass transfer alone and thus provides 
no immediate help in identifying the exact nature of these two sys- 
tems. However our hydrodynamical stabilization mechanism nev- 
ertheless plays a significant role in the orbital evolution of semi- 
detached ultra-compact binary white dwarfs. It must therefore be 
taken into account for identifying the region of parameter space 
where these systems are stable to mass transfer. 



1.4 Outline 

Our paper is structured as follows. We start in section|2]by briefly 
reviewing the essential material from perturbation theory of ro- 
tating stars required for our analysis. We give a more detailed 
review in Appendix |A| for the reader interested in a streamlined 
introduction t o the detailed, heavily mathematical formalism of 
ISchenkefg/.! ( I2OO2I) . In section [3] we derive the evolution equa- 
tion for the spin frequency of the accretor including tidal effects. 
This evolution equation is essentially a version of ^ where Jtidai 
and J2[i, ^] are given explicitly in terms of a mode expansion of 
the fluid displacement vector field. We also give the evolution equa- 
tion for the orbital frequency, needed to justify approximations used 
later in section|4l In section|4]we construct a simple dynamical sys- 
tem where we assume that a single mode is dynamically relevant 
for spin evolution. There we solve the equations of motion of this 
dynamical system both analytically and numerically, modeling the 
accretor as an incompressible, slowly rotating Maclaurin spheroid. 



^ Unless the binary is in accretion turn-on phase, where the mass trans- 
fer rate can be substantially below its equilibrium value, in which case it 
possible to have overall increasing orbital frequency. 



2 ELEMENTS OF PERTURBATION THEORY OF 
ROTATING STARS 

In this paper we use perturbation the ory of uniformly rotating stars 
as presented by ISchenk et al. which is very mathematical 

body of work. In this section we present the essential elements re- 
quired for the analysis of this paper, and give a more detailed re- 
view for the interested reader in Appendix IaI We use geometric 
units G = c = 1 throughout this paper. 

The perturbed stated of a uniformly rotating star can be en- 
tirely described by a Lagrangian fluid displacement vector field 
£,(x,t), which tells how fluid elements are displaced from their 
background position x in the frame corotating with the star. This 
fluid displacement is expanded in phase space as follows 



tt{x, t) 



■Ka{x) 



(25) 



where tt is the momentum canonically conjugate to ^. A given 
mode with associated eigenfrequency Ua is solution to a spe- 
cific eigenvalue equation [Eq.(|A9}]. The mode amplitude coeffi- 
cients Ca{t) obey the following evolution equation, which is de- 
duced by substituting expansion M5\ into Euler's equation. 



Ca+iLOaCa = -^ / (f X p $_l,{x) ■ ac^t{x ,t) , 
Oct 



(26) 



where p is background mass density and Cext is the acceleration 
field experienced by the fluid elements due to external perturba- 
tions, e.g. a companion object in a binary. The constant ha is given 
by 



(27) 



where Q is the angular frequency of the star and f2 is a unit vector 
pointing into the direction of the background star's angular mo- 
mentum vector. The complete tidal response of a rotating star to 
a prescribed external perturbation is then obtained by solving l |26l l 
for all modes of the star and substituting the results in expansion 
(|25}. 



2.1 Adiabatic approximation for time-dependent spin 
frequency 

So far the formalism discussed here applies for constant back- 
ground spin frequency ft. In accreting white dwarf binaries, the 
spin of the accretor will change in time so we need to discuss how 
to incorporate this in our analysis. A few complications arise in the 
case of time-varying spin. However if the spin of the star varies 
slowly in time compared to eigenfrequencies of dynamically rel- 
evant modes, we can still use normal modes that are solution of 
l |A9t in the following way. Consider the (infinite) sequence of vec- 
tor spaces of solutions to JA9l l for some range of spin frequen- 
cies (fimin, f^max), where riinin > and rimax < ribroak-up- 

Label the normal modes as {^^{x,n),LJa{Q,)), where Q denotes 
which vector space the mode labeled with quantum numbers a 
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belongs t(jf[ Assume now a displacement of the form ^{x,t) 
exp [ — j / ijJa(^)dtj ^ai^i We then have 



(28) 



the second term being a measure of the change in the mode eigen- 
function as the star is spun up or dowr0 . If the timescale of change 
of the spin frequency is much longer than the inverse normal mode 
frequency, the second term is negligible and we get 

,duJa 



2 t 



(29) 



Again if the timescale of change of the spin frequency is much 
longer than the inverse normal mode frequency, we can drop the 
second term in the right-hand side and we then see that our ansatz 
for the fluid displacement satisfies JA9I ) with B up to 

corrections of order Q/{Qujct), which we assume are much less 
than order unitjQ In this adiabatic-type approximation, we can 
then expand a generic fluid displacement in phase space as follows 



Tr{x, t) 



+ c.c. 



(30) 



where the set {£^^{x,Q.),iia{x,Q.)} are the normal modes and 
their conjugate momenta obtained from perturbation theory of stars 
spinning at constant frequency. The time evolution of the mode am- 
plitude coefficients is then given by ( |A25t , again up to corrections 
of order Q,/{yiuja), where now and ba depend slowly on time 
due to evolution of spin frequency. In the rest of this paper we will 
always omit the labels Q. on the modes and their frequencies. It will 
be understood that a given mode {S^^,uja) depends on spin fre- 
quency in the context of the adiabatic approximation detailed here. 



2.2 Mode amplitude evolution equation for tidal excitation in 
binary systems 

We now specialize to the case where the external acceleration vec- 
tor Oext is generated by the Newtonian gravitational field <I>ext of 
a companion object. We model the companion simply as a point 
particle of mass A/2 and assume that the orbital angular momen- 
tum is aligned with the spin axis of the perturbed star. The mode 
amplitude forcing term [the right-hand side of ( IA25H is then given 
by 



/ (fx spa $oxt. 



(31) 



where Spa is the Eulerian density perturbation of mode a. Since 
the background star is axisymmetric we may write the density per- 
turbation due to mode a as 



Spa = ga{r,6)e 



(32) 



In this paper we will use the convention > 0, which implies 
LOa may be positive or negativf^ 

Writing the orbital separation vector d pointing from the ac- 
cretor to the donor as d = d{cos (jiorb, sin ^orb, 0). We then make 
use o f the following ex pansion for the external Newtonian potential 
'I'oxt fao&Lail l fT99^ ) 



"™'^"-'"y!m(e,<^), (33) 



where x are inertial coordinates and where 



(! + m)/2 



47r 



{l + m)\a-m)\ 



1/2 



(34) 



the symbol (— )" being zero if n is not an integer. Substituting ex- 
pansion l[33) into l |31| l and using (f) = (p — J Q,{t')dt' then yields 



fa = Tae^P —irUa y(porb{t) — J Q{t')dt' 

where is the following complex number 



(35) 



ha 



The mode amplitude then satisfies 



(36) 



(37) 



If the force term oscillates at a frequency far from the mode's nor- 
mal frequency, we may write the following approximate solution to 



rrlaU — LOa 



(38) 



Equation J38b is only valid if the mode has not encountered a res- 
onance at an earlier time, in which case ([38} needs to be supple- 
mented by a homogeneous solution to ( 1321 obta ined by ma t ching to 
the solution to ([37} during resonance ((Flanagan & Racine! ( l2006h ). 



° One does not necessarily need to consider this sequence of Hilbert spaces 
to analyze the case of time-dependent spin frequency. We refer the reader to 
Appendix A of Bondarescu, Teukolskv & Wasserman (2007) for a detailed 
discussion of perturbation theory of stars with generic time-dependent spin 
frequency. 

^ It is possible to make a continuous, one-to-one identification of modes 
between vector spaces of neighboring frequencies, as long as one is a finite 
distance away from £7 = 0. The change in mode eigenfunctions can be 
meaningfully computed with the help of this identification. 

For compact white dwarf binaries with accretion rates of order 
10~®Mq yr~^, spin frequency of order orbital frequency and mode fre- 
quency of order spin frequency, the ratio tl/Q,u)a is of order 10^** — 10~^. 
The adiabatic approximation is thus fully justified for such systems. 



3 EVOLUTION OF SPIN AND ORBITAL FREQUENCIES 

In this section we give evolution equations for the spin frequency of 
the perturbed star and the orbital frequency of a binary undergoing 
conservative mass transfer. These quantities are essential in solving 
the mode amplitude equation of motion {37} since they determine 
the phase u{t) of the forcing term. 



The convention used in ISchenk et al\ j2002l) is to choose uJa positive 
and thus differs from ours. If needed, see AppendixlAlfor more details about 
mode pair notation. 



Tidal coupling in white dwarf binaries 7 



3.1 Evolution of spin frequency 

We obtain the evolution equation for the spin frequency from the 
fundamental equation 



dJsi 



dt 



(39) 



where T^^t is the external torque acting on the star. In Appendix 
|B]we show explicitly that the total angular momentum of the per- 
turbed star is given by 



Jstar — Iq^O + 



1 



din/, 



din Q.0 



lQm+ / d'^spfi- X 7r),(40) 



where Jo is the unperturbed star's moment of inertia, fio is the 
unperturbed star's spin frequency, and 50 is the perturbation in 
the star's spin frequency, so that the full spin frequency is = 
f2o + Assuming that the change in the star's moment of iner- 
tia due to changes in spin and mass (because of accretion) can be 
neglected, l|39) then gives 



dQ. 
HI 



Text - ^ / d^xpQ.-{£,x tt) 



(41) 



where / now stands for the moment of inertia of a non-spinning 
star, as we neglect the effect of spin on I. Equation | |4U is the de- 
sired evolution equation for the spin frequency. In Appendix IB] we 
derive a form of the spin evolution equation in terms of mode am- 
plitude coefficients, which is given by Eq.l[BT5) below. 

At this point of our analysis it is worth discussing the com- 
patibil ity of Ea. iB12t with the theorem of .Goldreich & NicholsonI 
l ll989h . who show that the secular change in specific angular mo- 
mentum of a fluid element following a suitably defined "back- 
ground trajectory" vanishes if the star is subject to the gravitational 
perturbation of a companion object rotating around the star. We 
want to point out here that this theorem does not necessarily imply 
that the spin of the star cannot be altered, in the sense discussed in 
section lA2l by tidal interactions. This important point can be illus- 
trated simply as follows. Label by xo the location at t = of a 
fluid element of the unperturbed star. We denote the trajectory this 
fluid element would follow if the star were to remain unperturbed 
at all times by x{xii,t). Assume now that at t = a perturba- 
tion is turned on and that the exact state of the star is described by 
fluid elements following the trajectory y{xQ,t). In the language of 
perturbation theory described in this paper, the fluid displacement 
vector is defined as 



y{xo, t) = x{xo, t) + I (a;o, t). 



(42) 



Note however that the traje ctory x(xo,t) does not descr ibe the 
background state defined in IGoldreich & NicholsonI fT989h . They 
instead decompose the exact state as follows 



y{xo,t) = x*{xo;t) + $,*{xQ,t), 



(43) 



where the time average of the perturbation ^* about the background 
trajectory x* is required to vanish. Now in general, the time aver- 
age of the perturbation ^ about the unperturbed trajectory x does 
not vanish, one specific example being when a mode of the star is 
resonantly excited. Hence what Goldreich and Nicholson call the 
background state in general differs from the unperturbed state. To 
linear order in ^, we have 



x''{xo,t) = x{xo,t) + {$,{xo,t)) + 0(r), 



(44) 



where here the brackets { ) denote a suita bly defined time averag- 
ing procedure about the unperturbed state ( IGoldreich & NicholsonI 



(fl98?)). The point is that the quantity (^) contains in general differ- 
ential rotation modes and therefore the background state described 
by fluid worldlines x* does not generically describe a fluid in uni- 
fo rm (rigid) rotation. One therefore may not interpret the theorem 
of lGoldreich & NicholsonI jl 9891) as implying that tidal interactions 
cannot change the s pin of a star, in the sense define d in section lA2l 
Still the theorem of lGoldreich & Nicholsoj d 19891) may in prmci- 
ple be converted into an evolution equation for the secular change 
of the spin frequency if one decomposes the background state onto 
the basis of Jordan chain modes describing arbitrary differential ro- 
tation and then isolates the rate of change of the coefficient ci of 
the uniform rotation Jordan chain mode. 

Now in the special case where {$) vanishes, then the back- 
ground state of Goldreich and Nicholson does correspond to the 
unperturbed state (which here is a uniformly rotating star). In that 
case their theorem does indeed imply that the secular change in the 
spin frequency of the star vanishes, in the context of a star perturbed 
by an orbiting companion object. We can easily show explicitly 
that this result is compatible with the equation of motion ( IB12b we 
derived rigorously using perturbation theory of uniformly rotating 
stars. To carry out this proof one first makes the ansatz f2 = to 
solve the mode amplitude evolution equation J37b . One then finds 
that CA{t) = \cA\er^"^^^^^\ where \ca\ is constant. Clearly the 
time average of CA(t) over a timescale long compared to l/ii is 
zero and therefore (4) = 0. Using this form of the mode amplitude 
coefficients, the form l IDlOb for the total external torque applied on 
the star and averaging Eq. JB12t over a timescale long compared to 
1/u, one finds immediately that the time averaged change in the 
spin frequency vanishes, consistent with the ansatz used to solve 
the mode amplitude evolution equation. Thus a constant spin fre- 
quency is the solution to the time averaged spin evolution equation, 
consistent with the theorem of Goldreich and Nicholson. However 
in more general situations where the frequency ii of the driving 
term in the mode amplitude evolution equation J37b varies in time, 
leading possibly to resonant excitations of normal modes, one ob- 
tains the time evolution of the spin frequency of the star by solving 
simultaneously l l37b and l lB12b . as long as the adiabatic approxima- 
tion of section[2T1is valid. 



3.2 Orbital frequency evolution 

We next give the evolution equation for the orbital frequency uIoa- 
Below we assume a Newtonian quasi-circular orbit, i.e. Wo^t, = 
(Ml + M2)ld^ = M/d^ with M being constant. We follow 
iMarsh. Nelemans"& Steeghsl ( l2004l) closely. We assume the orbital 
angular momentum Jarb evolves due to three effects, namely emis- 
sion of gravitational waves, mass transfer and tidal coupling. We 
may then write 



Jorb ~ Jgw + Jacc + Jtid 

where 

^GW — — ^^/iAf^ Jorb- 



(45) 



(46) 



In Appendix Id] we review the computation of Jtidai in Newtonian 
gravity and show that it can be written to leading order in the fluid 
perturbation as a sum over modes as follows 



d'tidai = — ^ mtjba^lcal^ + ~ i^oih) 



(47) 



the second term in the right-hand s ide modeling dissipa tive tidal 
coupling. As noted in Appendix K of lSchenk et a/.l ( l2002l) . each in- 



8 E. Racine, E. S. Phinney and P. Arras 



dividual term of the above sum does not correspond to the angular 
momentum deposited in mode a; the term J2[C,C] [cf- Eg.l lBlH 
giving the angular momentum in the star carried by the perturba- 
tion contains in general cross-terms between different modes. For 
a quasi-circular orbit and conservative mass transfer we have 



J O] 



M2 
M2 



(48) 



where q is the mass ratio K'h/Mi. We then obtain the following 
evolution equation for orbital frequency 



J or 



32 ,,2^Jacc ,^ .M2 



■ > ■maba — \Ca\ + 



dV 



JorbTS 



(49) 



Comparing l lB12b and l |49t and assuming that the binary is in a 
regime where Q. ~ Worb, i-e. close to synchronization, we then see 
that roughly, Coovh ~ {R/dftl. Thus for the purpose of solving 
( I37t . we will assume as a first approximation that uJoA is constant, 
since u{i) is dominated by spin frequency time evolution. We shall 
leave including back-reaction on orbital parameters when solving 
for the mode amplitude for a future publication. 



4 A SINGLE GENERALIZED ROSSBY MODE AS 
SYNCHRONIZING MECHANISM 

In this section we develop and analyze a simple dynamical system 
that attempts to extract the dominant features of the synchroniza- 
tion mechanism we propose, namely the stabilization of the accret- 
ing star's spin frequency through resonant tidal excitation of nor- 
mal modes. The essential features of the physical system we model 
in this section may be summarized as follows. To begin with con- 
sider a compact, detached binary white dwarf system. The binary 
is compact enough that it coalesces due to gravitational wave emis- 
sion. The system eventually reaches a point where the larger, less 
massive white dwarf fills its Roche lobe and mass transfer begins. 
The accreting white dwarf spins up as it accumulates material from 
its companion and eventually reaches a spin frequency at which 
one of its normal modes is resonantly driven. The amplitude of the 
mode then grows rapidly and may affect significantly the evolution 
of the spin frequency, depending on the size of mode back-reaction 
[cf. Ea.l lBT2t l. 

Our aim is to solve the coupled dynamical equations l |37t and 
( IB12t for mode amplitude and spin frequency for a single resonant 
mode and use this solution to catalog which modes are suitable can- 
didates for stabilizing the spin of the accretor. In this paper we will 
focus s olely on the excitation of generalized Rossby modes dBrvanl 
l ll889l) : iLindblom & Ipseil lm% \ since their normal frequencies 
scale as the star's spin frequency, and may then resonate when the 
star's spin frequency is of order orbital frequency. Other modes like 
/ and p modes resonate only at spin frequencies much larger than 
synchronized spin frequency. It is then likely that dissipative tidal 
coupling prevents the accretor to reach a large enough spin for such 
modes to stabilize it. For stars with buoyancy, g-modes are also po- 
tential candidates for spin stabilization at spin frequency of order 
orbital frequency, but we will not investigate them in this paper. 
In Appendix |C] we review the properties of generalized r-modes 
needed for our purposes. For a modem, in-depth discussion of these 
modes, the interested reader is invited to consult lEindblom & Ipsed 
l ll999l) . 



4.1 Main assumptions and evolution equations 

The dynamical system we study in this section is of course a crude 
approximation to the true dynamics of the binary. The main as- 
sumptions we make when solving l |37t and l lB12t are the following. 



(i) First we assume, as mentioned at the end of section 1X21 that 
orbital frequency is constant for the purpose of solving the mode 
amplitude evolution equation J37t . 

(ii) We assume the accretor spins slowly enough so that (a) we 
may use appropriate approximate formulae for the r-mode eigen- 
functions and eigenfrequencies and (b) we may consider its mo- 
ment of inertia as constant. 

(iii) We assume a constant torque on the star due to accretion of 
material from the companion. 

(iv) We assume that the Rossby modes are driven solely by the 
gravitational tidal field of the companion, i.e. we neglect the exter- 
nal forcing due to impact of infalling material on the surface star 
when computing the driving term in Eq.l|37|l. 

(v) Lastly we assume that the evolution of the spin frequency 
described by Ea.ll41t is dominated by a single r-mode. 

A given r-mode may be a suitable candidate for tidal synchro- 
nization only if it satisfies the following four criteria: (i) its az- 
imuthal quantum number m must be non-zero, otherwise the tidal 
driving term does oscillate in time [cf. Eq.(|37H and resonance is 
not possible; (ii) its quantum numbers must satisfy 



I + m — even, 



(50) 



otherwise the mode does not couple to the Newtonian tidal field [cf. 
Eqs.ll3l,l|36ll and lICTetiR (iii) for m > the mode's dimension- 
less frequency w = uj/2Q must satisfy the condition 



A = m + 2w > 



(51) 



for resonance to be at all possible, as can be seen from i ll It . and (iv) 
the mode's back-reaction parameter Va, encountered in Eq.(|6l and 
defined precisely by Eqs. l |52| l and l |57t below, must be positive and 
large enough (how large will be determined later in subsection l4.2t 
for the mode to be a candidate for tidal synchronization. As shown 
in table lEll in Appendix |E] there are five candidate modes, with 
I < 6, satisfying criteria (i)-(iii). Each candidate mode is marked 
by an asterisk in that table. We shall see below that modes will 
/ > 7 are not expected to satisfy criterion (iv) for generic binaries, 
which is why they are excluded from table lEll 

As mentioned above our model dynamical system consists of 
two coupled differential equations, one for the (complex) mode am- 
plitude Ca and one for the spin frequency of the star Q. For a single 
mode, spin frequency evolution equation ( I41t reduces to 



n = - 
1 

" 7 



Te,t - 2 (S7 X -iiuci„ + n X - 

Text 



d I |2 



d , .2 



(52) 



Conservation of angular momentum yields 



~ — Jtid 



Note that this condition excludes the so-called classical r-modes since 
they have Z = m + 1. This restriction on Z -|- m is lifted if we allow mis- 
alignment of the spin and the orbital angular momentum. However since we 
expect most of the white dwarf spin to come from accretion, the amount of 
misalignment should be small and tidal coupling of r-modes with I + m 
odd should be suppressed by some power of this small misalignment angle. 
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'^^a ba. I | J acc 5 



(53) 



where Jtidai and Jacc are the rates of change of orbital angular mo- 
mentum due tidal interactions and accretion respectively. We then 
obtain the following system of equations 



^ — ^acc i^a^^lCal ? 

where the quantities ilacc and Va are given by 



^acc — j^J acc ; 



(54) 
(55) 

(56) 
(57) 



where tl^cc is assumed constant. For the analysis of system i54i - 
i55\ performed in the next subsection, the only information we 
need about the parameters J^a and is that they are of the form 
(constant) x SI. Explicit expressions for JT^ and i?a for generalized 
r-modes of slowly rotating stars are provided in Appendix IeI 

4.2 Approximate analytic solution 

Let us start our analytic investigation of system (I54t-(I55I) by defin- 
ing a few quantities as follows 

7q 



Xa = ma + 2Wa, 
I la = -r — UJoih- 



(58) 

(59) 
(60) 
(61) 



Definition l |59| l is such that parameter ija is dimensionless and con- 
stant. We can then rewrite the system ll54b-ll55t as follows 



dt 



7Q + 



iXa{n - x[i + o{n/n^)] 



(62) 



(63) 



The scaling of the error terms in l l62t comes from assuming that the 
timescale for evolution of \ja\'^ is at most the "no back-reaction" 
resonance timescale trcs ^ Cl^^^'^ . The main advantage of this for- 
mulation of the evolution equations is that we have gotten rid of the 
term e"™"" in the mode amplitude equation of motion. However 
the presence of the Q"* factor in the evolution equation l |62t for the 
spin still makes the equation very complicated to solve. But since 
we expect back-reaction of the modes on the spin to be dynami- 
cally relevant mostly in the vicinity of the resonance, we simplify 
i62\ by replacing by Q]^^ in the second term of the right-hand 
side. We can estimate the errors introduced by this approximation 
as follows. Away from resonance the solution for the rescaled mode 
amplitude 7^ is given by 



1 + 



n 



+0 



yielding simply 

ti = n^,4i + o{ua)] 



(64) 



(65) 



in that regime. Numerically the 0{ua) corrections lie within the 
range 10~^ — 10~^ for modes of interest, independently of whether 
or not Q, is replaced by fij^"'' in the right-hand side of l |62| l. Thus re- 
placing Q. by Jljj'^^ in the second term of the right-hand side of l l62t 
introduces negligible corrections away from resonance. Now ap- 
proximations l |64| ( and ([65} fail when resonance regime is reached, 
which is when the dominant error terms in l l64t become of order 
unity, i.e. 



or alternatively when 



1 + 



(~)rcs 



(66) 



(67) 



For accreting binary white dwarfs these correction terms are typ- 
ically of order 10~*. Thus approximating Q, by QI^^^ in the right- 
hand side of l |62| ( introduces negligible corrections also in the reso- 
nance regime. Defining 



-^rcs ^ 3 



(68) 



the system of equations we solve here can be written in its final 
form as 

d , 

+ iXa [n — ^^a^\ ' 



7. 



7. 



= 1. 



(69) 



(70) 



The initial conditions we impose to integrate system J69l>-J70t are 
the following. We assume accretion begins at t = and that the 
accretor is not spinning prior to accreting phase. This implies of 
course f2(0) = 0. For the mode amplitude we assume that it starts 
at its equilibrium tide value, namely 7q(0) = ijXaS^a^ , obtained 
by setting 7a = at t = in J70t . We are now ready to solve 
(|69l)-l|70b analytically. Writing mode amplitude 7^ as 



la 



R , ■ I 

la + i7 



(71) 



with 7a real, we may rewrite l |69l l using l |70t as follows 

tl = fiacc - 2ua'yS- (72) 

Taking a time derivative of i72\ and using l |70| l once more yields 



Q + 2va ~ —ha (n — f^Iir^) , 
where the quantity ha is defined as 



ha — 2lyaXr 



:7q 



Then the imaginary part of l |70| l may be rewritten as follows 



ha 



n). 



(73) 



(74) 



(75) 



So far all that has been done is rewriting system l|54t-l|55ll in a 
slightly different, but very useful form. We are now ready to iden- 
tify the crucial approximation needed to solve l |73l l and f!5\ simul- 
taneously. The key lies in realizing that for a large enough value of 
i^a, which will be determined below, we may discard the fl term 
on the left-hand side of l |73l l. This approximation essentially yields 
some average spin frequency since we are dropping a term that 
generates oscillations in fl at frequency ~ \/Ti^. This average fre- 
quency solution is obtained directly from Eq.l|73ll and is equal to 



2Ua 

ha 



Substituting ( 176b into ( 175b then gives 

• 2X^Ua^a.cch^ 



hi +A\l,v'i 



(76) 



(77) 
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which is then easily integrated. The result is the following cubic 
equation for ha 



hi 



ha — ai{t)ha 



ao 



- 0. 



(78) 



The average frequency is then given by i76b . choosing the unique 
root of ( I78t that satisfies the initial condition ha{0) = 2h'a/i^a^- 
The explicit form of the complete solution is 



n{t) = fir, 



2j^a 



ao 



+ 



27 



1/3 



fflo 
2 



1/3^ 



27 



where the critical time tc is defined such that 



aiitc) = 



(79) 



(80) 



Clearly solution l l79t is not particularly revealing. In the large t 
limit however solution l |79t assumes the following very simple form 



where 



Aa^acc to) 



1/2 



for t » to, 



to = 



^acc 



(81) 



(82) 



The constant to is essentially the time it would take to spin the star 
up to the resonant frequency in the absence of mode back-reaction. 

It should be clear that approximation M6\ . leading to the large 
t solution J81t . must necessarily fail when back-reaction param- 
eter Ua is small, since in that case the star spins up through the 
resonance. An argument that shows this goes as follows. Assume 
( I76t and | |77I > to be valid aXt = Q. Initially we have /i^ > and 
hence ha > 0. Equation l l77t then implies ha > ha (0) > for 
all time, which means that Q < for all time, which cannot 
be true for Va less than some critical value Vait, since in the limit 
fa 0, the solution for the spin approaches — * Cl^cct con- 
tinuously. Approximation l |76t can therefore work well for all time 
only for modes that will stabilize the spin frequency, i.e. modes 
with i^a > J'crit. The value of i^crit may then be estimated by find- 
ing when approximation f!6\ fails, which is when fl ~ 2i'a- Using 
f!6\ and l l77t . we have 



n 

2v„ 



-3(2A„i^af^a 



h% 



{hi + 4a; 



2,y2X3 



(83) 



This expression has a maximum at ha = 16A^i'q/5. Putting this 
value for ha in | |83I ) and requiring f2 ~ 2ua then yields the follow- 
ing value for the critical value Vcrit of the back-reaction parameter 



I^crit = 0.24V Acf^ 



.3/2 
'acc ■ 



(84) 



Along with the late time solution ( 1811 ) for spin frequency, equation 
([84} is one of our key results. Again only modes with Va ^ fc-cit 



may effectively stabilize the spin of the accretor through tidal reso- 
nance. Modes with back-reaction parameters Va ^ I'crit do not have 
a strong enough stabilizing effect on the spin to prevent the accretor 
from being spun up through the resonance. For the candidate modes 
shown in table IE II %/A^^ 1.7 so a more conservative form of ( I84l l 

'3/2 

that we may apply to all candidate modes is i^crit = 0.4 Slacc ■ An- 
other useful version of bound ( |84t that must be satisfied by spin 
stabilizing modes is the following 



^ » 0.24v/a„ 



Ta 



i La 



-2 /61/2 

^acc 



i in 



(85) 



The latter version of the critical bound highlights better its depen- 
dence on tidal coupling strength, as Vaj^a^ is independent of bi- 
nary parameters while the back- reaction parameter Va, to which 
bound ( I84t refers to, depends on those parameters. 



4.3 Numerical results 

Here we present results of numerical integration of the system 
of equations ( I54t -(l55t and compare with the approximate ana- 
lytic solution derived in section [412] above. For realistic accretion 
rates, numerical integration of ( I54t -(l55t is challenging because 
of a drastic separation of timescales. The timescale of accretion- 
induced spin-up is of order ~ 10~^ yr~^5 while the characteristic 
timescale of the mode amplitude equation away from resonance is 
^ ~ 10~^ s~^. Stable numerical integration of system ( I54t - 
dSSi l using a Runge-Kutta scheme from S7 = to SI = fires takes a 
very long time and is not a desirable approach. What we do instead 
is start numerical integration close to resonance and use initial con- 
ditions for mode amplitude derived from the approximate solution 
( I38t to equation of motion ( I37t . In figure|2]we show the time evolu- 
tion of the spin frequency of the accretor obtained from integrating 
(I54il-(i55ll using such approximate initial conditions. The resonant 
mode used here is the Z = 5, m = 3, u) = —0.5266 generalized 
r-mode, and we used the same values for tidal ratio and mass ratio 
as reported in table IE II We normalize units of time fires = 1. In 
these units we used a value of fiacc ~ 10~^, its value in Hz/yr de- 
pending on the value of the orbital frequency. For an orbital period 
of ~ 5 minutes, the accretion spin-up rate is • 
For these parameters bound ( I85t is 



3 X 10"^ Hz/yr. 



;^>2.2xl0"^ (86) 

ores ^ ^ f 

^'53 

which is satisfied as one can see from table lEll We thus expect the 
mode to stabilize the spin for this particular choice of parameters. 

In figure [3] we show the relative error between the spin fre- 
quency obtained from numerical integration of J54t-(l55t and the 
full analytic approximation ( |79i l. The relative error is less than 
10~^, which gives additional strong justification to the various ap- 
proximations made in deriving this analytic solution. We attribute 
the initial discrepancy between the numerical solution and the an- 
alytic solution to our practical choice of initial conditions for nu- 
merical integration. 



4.4 Lifetime of synchronizer 

The synchronizing mechanism we propose is based on exciting a 
mode close to its resonance frequency over a long period of time. 
This process obviously cannot last forever since the mode ampli- 
tude grows as ^/t — to in the regime where the spin is stabilized, as 
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1.5 X 10 



Figure 2. This figure shows the time evolution of accretor's spin fre- 
quency obtained by integrating system )54t - j55t numerically. The follow- 
ing model parameters were used for moment of inertia of the accretor, the 
tidal ratio, the mass ratio and accretion-induced spin-up rate respectively: 
I = 0.4Mi?2, R/d = 0.2, q = 0.25 and f2acc = 10"^. 




2e+06 4e+06 6e+06 8e+06 1e+07 



2e+07 1.4e+07 



Figure 3. This figure shows the relative error between the spin frequency 
obtained numerically and the approximate analytic solution derived in sub- 
section |T2l 



can be seen by combining {STJ and ( 

hc\ = |7q(0)1 h 

^acc 



{t-to). (87) 



After a certain amount of time the mode amplitude will reach 
a value where the resonant mode couples significantly to other 
mod es of the star thro u gh non-linear hydrodynamical effects 
(see ISchenk et al\ j2002l) : iBrink. Teukolskv & WasserniM] j2004 



I2OO5 )), which saturates the growth of the mode amp litude. As 



^ mp li 

shown by ISrink. Teukolskv & Wassermanl ( l2004l [20051) , a precise 
understanding of the evolution of a perturbed star in the regime 
where a mode saturates requires a simulation involving thousands 
of modes. 

Denoting the saturation amplitude Csat , we obtain the follow- 
ing estimate for the duration Tm of the stabilization mechanism 
presented 



Tiif 



Csat 



T |20 

Q ^ ''a 



1 Csat 1 



10 I Csat I ^ 
^^rcs 



(88) 



for the / = 5, m = 3 r-mode and the same model parameters as 
used to generate figure[2] The lifetime of the linear synchronize^^ 
thus depends crucially on the mode saturation amplitude and could 
be as short as a few orbits but also as long as millions of years. 
For example an optimistic saturation amplitude of ~ 10^^ and a 
slow accretion rate, say 1O~^M0 /yr yields a lifetime of order 10® 
years. For sake of comparison, the radiation reaction timescale of a 
binary with components of 0.6 Mq and 0.2 Mq and orbital period 
of 10 minutes is frr ~ 8 x 10^ years. After the resonant mode 
saturates, the fate of the binary remains an open problem. 

4.5 Polytropic stellar model 

Clearly the incompressible model we used above is not very realis- 
tic. It is therefore worthwhile checking if changes in mode eigen- 
functions due to using a more realistic stellar model affect signifi- 
cantly the strength of tidal coupling and spin back-reaction parame- 
ters of our candidate synchronizers. Here we report values for mode 
frequencies, tidal coupling and back-reaction parameters computed 
from a polytropic star of index n = 3/2 (or equivalently Fi = 5/3) 
and spinning at frequency = Q.QIM / B? . There is no buoyancy 
in the star, so that g-modes are degenerate, zero-frequency modes. 
For this stellar model, the mode frequencies, tidal coupling and 
back-reaction parameters are, for the two main candidate modes 
discussed previously. 



(TOi=4,m 


= 2)n=3/2 


(Wi=4,m= 


-2) incomp 


(W; = 5,m 


=3)n=3/2 


(Wi = 5,m = 


-3) incomp 


\J~l^i,m-- 


=2|n=3/2 




-2 incomp 


\'^l — ^,m- 


= 3|n=3/2 




-3 1 incomp 


and 






= 2)n=3/2 


(t';=4,m= 


2 ) incomp 




= 3)71=3/2 




3 )incomp 


For the I 


= 5, m 



= 0.861, 



0.823, 



0.0897, 



= 0.0700 



1.94, 



= 1.65. 



(89) 



(90) 



(91) 



(92) 



(93) 



(94) 



pling strength [cf. Eq.l|92H is such that criterion l |85t is not cleanly 
satisfied anymore. For an accretion rate fiacc = 10~^ (in units of 



rate of order 10 



'Mq yr ^, bound l |85t becomes 



> 4.7, 



'53 



From l |94t and table lEll we have 



Ores 
''53 



= 24.02, 



(95) 



(96) 



fia 



By "linear synchronizer" we here mean the regime where the mode am- 
plitude is small enough so that it can be evolved using the linear, leading- 
order mode amplitude evolution equation, neglecting non-linear mode- 
mode couplings. However the coupling between the tides and the uniform 
(rigid) rotation mode is still non-linear (quadratic in the mode amplitude). 
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which is only a factor of 5 larger than the critical value given in 
( |95t . Therefore one cannot say whether or not this mode will act as 
a synchronizer generically, since slight changes in accretion rate or 
tidal ratio may push the mode on one side or another of the bound. 
For example tidal coupling strength squared for an / = 5 mode 
scales with tidal ratio as (Ri/dY^. Thus a slight change in tidal 
ratio may easily increase or decrease the critical value in bound 
l |95t by an order of magnitude. 

Lastly the lifetime of the synchronizer is still given by the fol- 
lowing expression 



Tm ~ I Csat 



(97) 



so for the same accretion rate, the lifetime increases only by a factor 
of 1.65 [cf. Eq.J94|l1 compared to the incompressible case. 



5 CONCLUSION AND FUTURE WORK 

In this paper we suggest a non-dissipative mechanism that can sta- 
bilize the spin of the accretor in an ultra-compact binary white 
dwarf system. The effect that stabilizes the spin of the accreting 
white dwarf is the back-reaction of a resonantly driven generalized 
r-mode on the uniform rotation mode. For the model we analyze in 
detail we assume that only a single mode is dynamically relevant 
(see section |4] for more details). We also assume slow rotation of 
the accretor, constant accretion torque, constant orbital parameters 
and a perfect incompressible fluid equation of state. By integrating 
the equations of motion of our model both numerically and ana- 
lytically we show that pure hydrodynamics may stabilize the spin 
of the accretor. No dissipative effects are required for synchroniza- 
tion. However the lifetime of the synchronizing mechanism is lim- 
ited since it requires driving the mode continuously very close to 
resonance. Eventually the mode enters a non-linear regime where 
it couples to other modes, which saturates the growth of the res- 
onant mode's amplitude. At that point the synchronizing mecha- 
nism most likely shuts off. A correct computation of the saturation 
amplitude involves simulating a network of coupled modes follow- 
ing the footsteps of iBrink. Teukolskv & WassermanI ( |2004| . l2005h . 
a challenging task that we leave here as an open problem. More 
future work includes solving the full equations of motion, where 
the accretion rate and orbital parameters are computed and evolved 
rather than assumed to be constant external input quantities. In ad- 
dition if mode damping timescales turn out to be comparable to 
the lifetime of the synchronizer, then one should include damp- 
ing when solving the mode amplitude evolution equation. It would 
also be interesting to investigate the properties of the stabilization 
mechanism if the accretor experiences core crystallization, as the 
properties of the r-modes (i.e. frequency spectrum, strength of tidal 
coupling, etc..) existing in the fluid shell surrounding the core dif- 
fer from the case of a fluid spheroid. 
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APPENDIX A: PERTURBATIONS OF UNIFORMLY 
ROTATING STARS 

In this appendix we review basic materi al from perturbation theory 
of uniformly rotating stars developed in 'Schenk ef a/.] ( l2002h that is 
required for the analysis of this paper. This appendix is rather tech- 
nical but the formalism reviewed here is essential in correctly ana- 
lyzing the tidal response of a rotating star. We hope to provide the 
interested reader an accessible introduction to the b eautiful, albeit 
mathematically heavy, work o f I Schenk et al\ ( l2002h . We refer the 
reader to that paper for complete details. We employ units where 
G — c — 1 throughout. 

Begin by considering an unperturbed star of mass Af and ra- 
dius R uniformly rotating at constant spin frequency fl. The radius 
of the spinning star is defined as the radius of the sphere enclosing 
the same total volume as the star. The hydrodynamical equations of 
motion in the frame corotating with the star are 

dp 
dt 

Dt p 



+ V ■ {pu) = 0, 



(Al) 



(A2) 



where p is mass density, p is pressure, u is the fluid velocity, $ 
the self-gravitational potential and Ooxt is an acceleration field pro- 
duced by some external perturbation. The operator D / Dt is the 
usual convective time derivative. For the unperturbed star, p does 
not depend on time and the velocity u vanishes. The corotating 
frame is related to the inertial frame by the simple coordinate trans- 
formation (f> — ip ~ J fldt, where (f) is the corotating azimuthal 
coordinate and is the inertial azimuthal coordinate. If a fluid 
element of the unperturbed star located at x is displaced by the 
vector ^{x,t), the linearized hydrodynamical equations for the La- 
grangian displacement vector { take the form 



(A3) 



where the action of the operator i? on a vector is defined as 

B ■ { = 2f2 X {. (A4) 

For stars without buoyancy, the action of the operator C on a vector 
is given by 



^Sp + S<i> 
p^ 



-Sp + 



p 



(A5) 



where F is the adiabatic index of the perturbation, which we here 
take to be the same as the background star. The quantities Sp and 
5p are the Eulerian perturbations in the pressure and density respec- 
tively. In terms of the displacement vector the density perturbation 
is 

5p=-V-(pO- (A6) 

The gravitational potential perturbation 5$ obeys 

V^5<1> = 4tt5p. (A7) 

The fluid displacement equation of motion l lA3t can also be ob- 
tained from the following Lagrangian density 

Z: = ii ■ i + if ■ B ■ { - i{ ■ C ■ C + a,,t ■ ^. (A8) 

The normal modes of the star are obtained by assuming displace- 
ments of the form ^{x,t) = e~'"'4(a;) and no external forces. 
One must then solve the following eigenvalue equation for the nor- 
mal modes 
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iujB + C 



■4(x)=0. 



(A9) 



A mode {^^,uja) is solution to Eq. JA9b . Later in the paper we will 
use the notation wa = ujaII^ for dimensionless mode frequen- 
cies. Note that given a mode (^^, cja) with uja 7^ 0, there exists 
another distinct mode (^^^ijJb) solution to ( |A9t with = 
and UB = —UJA- Thus modes with non-zero frequency always oc- 
cur in pairs. 



Substituting ( IA16b -l lA17l into l lA14b and l lA15t yields 

■KA = -iuJA^A + ^B-i^, (A20) 
fTA = il^ATA - ■ TA- (A21) 

If we now assume that all mode frequencies are real, it is then 
possible to show the following ( this proof is a little involved; we 
refer the reader to Appendix A of lSchenk et al\ j2002l) for details) 



Al Phase space expansion formalism 

The eigenvalue equation ( |A9t satisfied by a given mode is peculiar 
in the sense that the operator acting on the perturbation depends on 
the eigenvalue uo. This property leads to the fact that although it is 
possible to find a complete set of solutions to ( |A9l l so that a generic 
fluid displacement may be expanded as 



ta 



(AlO) 



the evolution equations for the mode amplitude coefficients Qa {t) 
are not in general uncoupled from one another. To get around this 
problem one needs t o look at the fluid pertur bation in phase space 
(ISchenk et g/l j2002h ; lDvson & Schutj ( 1 19791) ). A general fluid dis- 
placement vector is expanded in phase space as follows 



C{x,t) 



where 



i{x,t) 
Tr{x, t) 



A 



iAix) 

■KAix) 



2^^ 



(All) 



(A 12) 



is the momentum conjugate to { obtained from the Lagrangian den- 
sity l lASt . Using iA3l the evolution equation for the phase space 
vector ^ is easily seen to be 



-IB 



1 

4B 





flcxt 



T-C+F.(A13) 



Thus the eigenvalue equation equivalent to l |A9l l in phase space is 



[T + iuoA] Ca = 0. 



(A14) 



An important point to note here is that operator T is not Hermi- 
tian, which leads to the existence of left eigenvectors Xa distinct 
from the right eigenvectors (which incidentally do not form a 
complete basis; to form a complete basis, one needs to inclu de all 
Jordan chain modes [see Appendix A of iSchenk et aZli l2002l) fora 
complete discussion]) that satisfy 



[t^ - iLj*A] Xa = 0. 

We will write the right and left eigenmodes as follows 



Ca = 



Xa = 



^A 

tta 

TA 
TA 



(A 15) 



(A16) 



(A17) 



Now it is always possible to choose Xa '^^^^ t° Ca '■^6 

sense that 

(XaXb) = (o-a,^b) + {ta,ttb) = Sab, (A18) 
where the inner product of two vector fields is defined as 



= J d^'xpix)e{x)-i'ix). 



(A19) 



bA 



(A22) 



where bA is a (real) constant. This constant is determined from 
( IaTsI ). which may be rewritten as 



bASAB = {4a>«-B • Is) + {uJA+l^B){$A,iB)- 



(A23) 



The orthonormality relation l lAlSt allows us to invert expansion 
JAIU to get 



CA = -i^Ai. + is ■ I + Tt) 



Finally, i[AT8) may also be used to extract from l |A13l l the following 
evolution equation for mode amplitude 



CA+iidACA = — (^Ai'^cxt). 
OA 



(A25) 



To summarize, the tidal response of a rotating star is described by 
a Lagrangian fluid displacement vector field ^{x,t) and its canoni- 
cally conjugate momentum tt (a; , t). These two vector fields may be 
decomposed in phase space into normal modes according to l lAl It . 
a given normal mode {^j^,u)a) satisfying the eigenvalue equation 
l |A9t . The time evolution of a given mode amplitude is governed by 
Eq.( |A25l l. 



A2 An important set of zero-frequency Jordan-chain modes 

As mentioned in the previous section the modes analyzed above 
do not form a complete basis on which one may expand an arbitrary 
fluid perturbation. To form a complete basis one needs to add all 
Jordan-chain modes. A Jordan chain of length pA associated with 
a given eigenfrequency uja is a set of vectors Ca o-' — — Pa 
satisfying 



[T + i^A]CA,a^CA,a^l, 

for 1 < cr < Pa and 
[T + iiUA] Ca,o = 0. 
The associated left eigenvectors satisfy 

[T^ -iuj\] XA,a = Xa,<t-1> 

for 1 < cr < Pa and 
[Tt - iuj\] Xa,o = 0. 



(A26) 



(A27) 



(A28) 



(A29) 



The gene ralized orthonormal ity condition for Jordan chains 
adopted in lSchenk et al\ ( l2002h is 



(XA,<T>Cs,r) = SABSp^,a 



(A30) 



The equations of motion for the associated mode amplitudes in the 
absence of external forcing can be shown to be (see Appendix A of 
ISchenkef a/1j2002h ) 



CA.CT + iiUACA,a = CA, <t+i, for a < pA 

CA,pA + ii^ACA,pA = 0. 



(A31) 
(A32) 
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The main reason motivating us to discuss Jordan chains in this 
paper is that all differential rotation modes must be Jordan chains 
of length one. One proves this as follows. The fluid velocity pertur- 
bation corresponding to a differential rotation mode is of the form 



/?=(«o,«o) + (ii,C-4i) 



(A42) 



di{x,t) 
dt 



Su{x) — r±5fl{r±)e,f 



(A33) 



where r± — r sin 9. The velocity perturbation being time- 
independent first proves that the mode is a zero-frequency pertur- 
bation. Now any normal mode is a Jordan chain of length p > 0, 
so we have (dropping A subscripts) 



5u{x) = ^c^(t)$^{x) 



(A34) 



The mode amplitude equations of motion l lA3U -l lA32t show that 
Ccr{t), in the case of a zero-frequency chain, is a polynomial in t 
of order p — a. Since 5u is independent of time but non- vanishing, 
then the chain must be of length one. In particular, the uniform 
rotation mode (5fl{r±) — 5fl — constant) which moves the star 
from a given background spin frequency f2 to a new background 
spin frequency + 5Q falls into this category. From Eq.( IA34t one 
immediately finds that for a generic differential rotation mode, the 
mode function is 



(A35) 



where 7 is a constant fixed by an overall mode normalization con- 
dition. The defining equations of Jordan chains ( IA26I ) and ( |A27t 
are then used to show that d is solution to 

C-^i = -B-^o- (A36) 

In the case of uniform rotation, we use the convention 

io = nxx. (A37) 

The complete expression for the uniform rotation Jordan chain is 
therefore 



Co 



Ci 



Ci 



(A38) 



(A39) 



with {(, and ^ j given by l lA37t and l lA36t respectively. Through- 
out this paper, when we say the spin of the star changes, we mean 
that the amplitude of the above uniform rotation Jordan chain mode 
changes. The amplitude of the uniform rotation mode defines the 
natural coordinate system corotating with the star in which stel- 
lar perturbation theory is applied. In the specific application we 
consider in this paper, it is the time dependence of the coordinate 
system corotating with the star due to the evolution of the uniform 
rotation mode that is responsible for the resonant excitation of nor- 
mal modes in an accreting white dwarf part of a binary system. 

Lastly, it will be useful later on to know the left Jordan chain 
associated with the uniform rotation mode. The defining equations 
( lA28t and l |A29t , along with l lA30t and l lA38] l-l |A39t give 



Xo- p 



Xi = -3 



where 



So 



(A40) 



(A41) 



The constant l3 ensures that the left and right eigenvectors given in 
jA40t -l lA4n i and jA38t -l |A39l l satisfy the orthonormality condition 
l |A30t . 



A3 Mode pair index notation 

As mentioned previously, modes with non-zero frequencies always 
occur in pairs. We may then rewrite the phase space expansion in- 
dex A as the pair of indices (a,e), e = ±, with the following 
identifications 



(A43) 
(A44) 
(A45) 
(A46) 

(A47) 
(A48) 



Equations ( IA43t to ( IA46t and JA23t then imply 

ba,+ = be 
be- = —ba- 



Since the fluid displacement and its conjugate momentum are real 
vector fields, Eqs.(|A43) to ( |A46t , ( lA24l l and l |A47l l-l |A48t yield 

Ca,+ = Ca (A49) 
Ca~ = c* (A50) 

If the perturbation is dominated by modes with non-zero frequency 
we may rewrite expansion JAl It as follows 

CA(i) 



ax,t) 

'7T{x,t) 



A 



$.4(33) 
TTAix) 



TTa{x) 



where c.c. stands for complex conjugate 



+ c.c. 



(A51) 



APPENDIX B: DETAILED DERIVATION OF SPIN 
FREQUENCY EVOLUTION EQUATION 

The differential equation governing the time evolution of the 
spin freq uency ft of the star can be obtained using results of 
ISchenk et al. l(|2002l), who rigorously derive the evolution equation 
of the uniform rotation mode amplitude to second order in stel- 
lar perturbation theory. Here we shall instead derive the evolution 
equation for the spin frequency following the procedure sketched 
in the in troduction. Our evol ution equation of course matches the 
result of ISchenk efoll 1 I2OO2I) . However we believe our derivation 
is slightly simpler to follow so it is worth including it here. We start 
from the expression for the component of the angular momentum 
of the perturbed star along its unperturbed spin axis, and expand it 
to second order in the perturbation. A direct computation yields 



J d''xp{x)n- ^^(x + i) 



Ho X {x + ^)+£, 



flo X x + n 
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+ 



j3 

a xp 



(flo X • (n X a;) + fl ■ (a; X tt) 



J^o/o + J2[C,C] + Ji[C], 



(Bl) 



where fJo — fioJ^ is the unperturbed spin frequency vector, which 
was denoted by f2 in Appendix IaI The full spin frequency of the 
star is now written as 



(B2) 



Let us now focus on the term Ji [C,] . Using the defining equations 
( IA37b and ( IA40t for the uniform rotation Jordan chain, we may 
rewrite Ji [C,] as follows 



j3 

a X p 



dPx p 



(f2o X 4) ■ (f2 X x) + n ■ (x X tt) 



■ (flo X lo) + 77 ■ ^0 



f3ci{t), 



(B3) 



where the last line has been obtained using Jordan-chain orthonor- 
mality condition l lA30t . Thus the piece of the angular momentum 
linear in the perturbation depends only on the amplitude ci of the 
uniform rotation mode. It is independent of all other modes of the 
star from orthonormality condition JA30t . It remains to relate the 
coefficient ci to the spin frequency perturbation and to compute the 
quantity f3 defined in ( IA42b . Using l lA38b . the velocity perturbation 



Su = r±S^l = 50, {ft x x) — SQ {q, 
which gives simply [cf. Ea. JA34b 1 

Co = ci = 5Q. 



(B4) 



(B5) 



We now turn to the computation of /3 [cf. Eg. JA42b 1. The first term 
of the right-hand side of JA42b is easily shown to be 



(B6) 



the unperturbed moment of inertia of the star about its spin axis. 
This result leads us to expect that the second term of l lA42b is re- 
lated to the perturbation of the star's moment of inertia under uni- 
form spin-up. The moment of inertia of the perturbed star about its 
spin axis is given by 



(B7) 



/ = J d\p {{x + if-[{x+i)-^f}. 

Expanding this expression to leading order in ^ yields after some 
algebra and an integration by parts 



d^xv ■ (p|)[a;^ - {x ■ ny 



(B8) 



If the perturbation consists of uniform spin-up, we then have, since 
p^o is divergenceless, 

/ = Jo - ci / d'xV ■ {p$.i)[x^ - {x ■ hf] = Io + 51. (B9) 



On the other hand we have, from l |A36b 

= ~ J '^^^ P^i ' [2^0 X (n X x)] 



flo J d^x ■ V [o;^ - (f2 ■ o;)^] 



Oo 



d^xV 



(p^ J [x' - (f2 . 



51, 



Cl 



where the last line has been obtained using 



Jo + ^51 

Cl 

d In Jo 



Cl 



1 + 



Iq50. 



(BIO) 
We finally have 



(BU) 



d In fio - 

In this paper we will assume that the star spins slowly enough so 
that the change in moment of inertia due to change in spin may 
be neglected. Taking a time derivative of jBlb then yields, in the 
context of the adiabatic approximation detailed in section l2n 



dO. 
HI 



Text - ^ / d^xpn - X tt) 



(B12) 



where Txt is the total external torque applied on the star, and / is 
now the moment of inertia of the non-rotating star. In the applica- 
tion we consider below there will be two contributions to the total 
external torque, namely an accretion-induced torque and the torque 
on the star due to tidal interactions. In terms of sums over modes, 
Eq.( IB12b can be rewritten as 



Text 

~r 

Txt 



1 d 

idi^' 

A,B 



C*aCb (rj X 



ILJB 



1 d 



I dt^ 

A,B 



CaCbOK-AB, 



(B13) 



where dAI lb has been used. Since oc e*™'*'^, the only terms 
contributing to the above double sum have rriA ~ niB, i.e. ICab oc 
5mAmB ■ modes that do not resonate, we may use approximate 
solution B8b for the mode amplitudes to split the double sum of 
( IB13b in two pieces, one that includes terms where either mode 
A or mode B resonates and another involving modes that do not 
resonate. We write this split as follows 



n = 



Txt 1 d 



I dt^ 

A.B 



1 ds-^ 
Idt^ 



(nr) 



caCbOKab 



Q.T*aTb 



{niAU — u;A){mAU — lub) 



ICab. (B14) 



Since all terms in the sum over non-resonant modes depend on time 
only through the star's spin frequencjl^, we may replace the oper- 
ator d/dt acting on each term of this sum by AabO/0, with A as 
being of order unity. This leads to the following expression 



n = 



(B15) 



where the effective moment of inertia I^b is given by 

l^{nr) /^abT'aTb 



1{\ 



A.B 



{vriAU — ujA){rnBU — ojb) 



Orbital frequency also depends on time but in this paper we will keep 
orbital parameters fixed for the purpose of computing the time evolution of 
the spin frequency. We motivate this approximation in the next subsection. 
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(B16) 

Although we are unable to prove that the sum over all non- 
resonant modes is small compared to unity, we suspect it is the 
case since a typical driving amplitude Ta scales as R/d lo some 
power p > 2, where R is the accretor's radius and d orbital separa- 
tion. Even though the binary is quite compact (i.e. R/d is not that 
small, maybe of order ~ 0.2), the leading term in the sum over non- 
resonant modes, which scales as (R/d)'^, is of order 10~^ — 10~^. 
It thus appears unlikely to us that the sum over all non-resonant 
modes may contribute to more than a few percent to the effective 
moment of inertia Jeff defined in ( IB16b . 



APPENDIX C: BRIEF OVERVIEW OF GENERALIZED 
R-MODES 

Here we briefly summarize the properties of generalized r-modes 
needed for this paper. This appendix re l ies he a vily on the work of 
iBrink Teukolskv & Wassermad j2004 l2005h : iLindblom & Ipsej 
( Il999h . Note however t hat we use the frequency sign convention 
of ISche nk et al.\ j2002|) instead of the frequ e ncy convention of 
[Brink, Te ukolskv & WassermanI j2004 l2005h : iLindblom & Ipsej 
( 119990 . For I > 2 the Lagrangian displacement of a given mode 
^A can be written as 



(CI) 



where tpA is a normalization constant fixed below and where the 
hydrodynamic potential 5 (7a is given by 



5Ua = Pa 



(0 Dm / - 



(C2) 



where I3a is a constant and where p., (j)) are a bispheroidal coor- 
dinate system defined below. We will henceforth drop all A sub- 
scripts for notational convenience since no ambiguity occurs in 
what follows. In cartesian coordinates, the tensor Q'-' is given by 



1 



w 



(C3) 



where w = uj/2Q and e'^^ = S'^'S^^ - S'^S^"". The bispheroidal 
coordinates p., tj)), esse ntial to write a sepa r ated s olution for SU, 
are explained in detail in ILindblom & Ipsed ( Il999l) . They are re- 
lated to cartesian coordinates by the following relations 



where 



[(1 + Co) 



(C4) 
(C5) 

(C6) 



(C7) 



the parameters a and i^o being the standard parameters describing 
the background spheroid, namely a is the focal radius and (o is 
related to the spheroid's eccentricity by = l/(l-|-(^o). Note also 
that bispheroidal coordinates depend on normal mode frequency. 
The volume of the star is covered by < C 1' ~Co < A < Co> 
where 



62 1 



(C8) 



We use the convention of iBrink. Teukolskv & WassermanI ( l2004h 
for mode normalization, namely a mode of unit amplitude has (ro- 
tating frame) energy equal to the fixed value MR'^Q,^. In the slow 
rotation approximation one can show (.Brink (.2005i) ) that the nor- 
malization constant ip is 



2Mi Q.^ 



1/2 



(/-m)! {21 + 1) 
(Z + m)! + 



1/2 



(C9) 



where Q. — f2/f2brcak-up = ^/\/ Mi/Rl- Lastly we give a for- 
mula for the density perturbation of a normalized mode 



Sp = i>p S{p) SU + 5<i> 



(CIO) 



where 5{p) is a delta function with the background pressure as its 
argument and 5$ is the gravitational potential perturbation of the 
mode. The density perturbation is non-zero only on the surface E 
of the star. On this surface the function 5$ is equal to 



imcj) 



(CU) 



w here -D(Co) is a kn o wn fu nction of i^o [cf. Eqs.(7.1) and (7.4) 
of ILindblom & Ipsed ( 1 19991) 1. It is also important to note that 
on the surface of the s tar, we have SU\j: = /3P,'"(/i)e""'^ 
jLindblom & Ipsed ( Il999l) ). In the slowly rotating limit, one can 
show 



2(/-l) 



(CI 2) 



The coordinate p appearing first in dCl lb is part of the spheroidal 
coordinate system defined by the relations 



z — a(^fi. 



(C13) 
(C14) 
(C15) 



In the slowly rotating limit, the spheroidal coordinate system be- 
comes the usual spherical coordinate system (r, 9,(j))on the surface 
of the star. In this approximation we then have 

21 + 1 



Sp = tpfJ 



2(1 ~ 1, 

M Q:^21 + 1 
R?~ l~ 1 



p'S{p)Pr{cos9) 



3(1 



1/2 



S{r - R)Yir^ 



7rZ(/ + l) 

(CI 6) 

where we made use of iC9h iClH and the slow rotation limit of the 
pressure distribution of an incompressible fluid spheroid, which is 

2np'^ 



P = 



2\ 

r ). 



(C17) 



APPENDIX D: TIDAL COUPLING TORQUE 

Here we derive the rate of change of orbital angular momentum 
induced by the excitation of normal modes in one of the stars of the 
binary. This result is well known but we include the derivation here 
in the context of perturbations of rotating stars for completeness. 
Consider a rotating star perturbed by the following Lagragian fluid 
displacement vector field 



${x,t)^"^CAmAix) 



(Dl) 



The net force Fcxt on the perturbed star due to the compan- 
ion's tidal field is given by 
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dMac^t{y,t), 



(D2) 



where dM is the mass of a fluid element and 1/ = a; + ^ is the lo- 
cation of that fluid element. Now we have dM — d?x p{x), where 
p{x) is the background density and x is the background location 
of the fluid element. To leading order in the perturbation we then 
get 



d^x p{x)aaxt{x + i,t) 



d X p{x) {ac^t{x,t) +4 ■ Vacxt(a;,t)} 



(1) f 



(D3) 



The zeroth order force term describes the binary's motion when no 
modes are excited. The term depending on ^ can be rewritten as 
follows 



-f'cxtIC] = J d'^x p{x)£_i{x,t)dia^^t{x,t) 

= ^ / (f'x p(x) CA(t)£,A{x)diaa^t{x,t) 



A 
A 



CA(t) I <fx5pA{x)ac^t{x,t). 



(D4) 



The external acceleration is given by 

aB^t{x,t) = —V^c^t{\x — d{t)\ 



\x — d\ 



(D5) 



where d is the binary's orbital separation vector and V'"*' denotes 
the gradient operator with respect to d, keeping the field point x 
fixed. We then have 



A 



(D6) 



where Olt has been used. If we write d = d(cos (p, sin (j>), the bi- 
nary's equation of motion (ignoring gravitational radiation reaction 
and mass transfer) is then 



d- d(ji -i — 



- icAbAddfX 
11. ^ — ' 



A 

-y 

ad ^ 

A 



iCAbAdrfifA 



(D7) 
(D8) 



The rate of change of orbital angular momentum due to tidal cou- 
pling is then given by 



>Aid 



A 



icAbAd4,fA 



A 



mACAbAfX 



A 



TUACAbA 



Ca — IUJACa 



(D9) 



where l |35l l has been used. Using mode pair notation, i.e. Eqs.l |A43"t 
to l |A46t . jA47t -l |A48l l and l |A49t -l |A50t , we may further simplify 
this result to finally arrive at 



d_. 

'dt' 



(DIO) 



APPENDIX E: COMPUTATION OF TIDAL COUPLING 
AND BACK-REACTION PARAMETERS 

In this appendix we gather analytic results for integrals required in 
section]?] 



El Mode amplitude driving term 

The first quantity we compute is the driving term J^a appearing in 
Eq.ll37t. To compute this quantity we need 



(i^.-V'I'ext) 



d X Sp'^^c 



(El) 



Using J33t and l lC16t we immediately obtain 



d 

3(1 



'4(r 



TVI{1 + 1) 



1/2 



(E2) 



The normalization convention we choose in Appen dix [C| implies 

ba = AIiRiQ/{2wa) (see also Appendix A of ISchenk et al\ 

| |2002|) ) and we get 



iWa 0,M2 



d 

d 



R^^'*\^^ (2^ + 1) /3(i-^a) 



2(/-l) V + 



(2/ + 1) 3(1-^1) 
^""2(7^V + ' ^^^^ 



where q — M2/M1 is the mass ratio. 

E2 Integrals in the spin frequency evolution equation 

We next compute some of the integrals appearing in Ea.l lB13t . 
These are more tricky since they cannot be reduced to contribu- 
tions from the surface of the star alone. It possible to evaluate these 
integrals analytically for a fluid displacement vector composed of 
a superposition of generalized r-modes using a recursion technique 
(Brink (2005)), but a general closed-form expression for the inte- 
grals appearing in Ea.l lB13t is not available. In the model where a 
single r-mode pair is assumed to be dynamically relevant however, 
the required integral is 

ifc,, = -iuj^iz X $^,$^) + il{z x$^,zx (E4) 

which we now give in the slow rotation approximation. Using l ICll l 
and iC3\ we have, dropping the a and 1 subscripts below for 
brevity, 



1604(1-^2)2 

X (^Xjk — 2iw ejkz 
xdp5U*dqSU, 



d\ S'" 



w' w 



(E5) 



where X 



Ojk ■ 



ZjZk- Performing the index contractions gives 



rp 



16^4(1 - 'U;2)2 



2 
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Table El. Numerical values for the parameter A = m + 2w, the resonant 
spin frequency r2ros, the tidal coupling \T\ and the back-reaction pai'ameter 



V for generalized r 


modes with / 


i- m even, 


up to Z = 6. 




{I, m, w) 


A n 


res /, , 




p/n 


3, 1, —0.755 


-0.510 




5.6/19e-05 


1.8537 


6, 1, U.Uooo 


1.177 


0.850 


1 .U076e-05 


-25.052 


4, 2, — U.olD 


0.768 


2.604 


5.9429e-06 


7.5596 


4,2,0.116 


2.232 


0.896 


1.4105e-06 


-38.004 


5,1,-0.903 


-0.806 




5.5717e-07 


0.84994 


5,1,-0.5228 


-0.046 




6.4004e-07 


5.1797 


5,1,0.0341 


1.068 


0.936 


4.8944e-08 


-62.488 


5,1,0.5917 


2.183 


0.458 


6.8504e-07 


-4.2286 


5,3,-0.5266* 


1.947 


1.541 


6.9443e-07 


14.566 


5,3,0.1266 


3.253 


0.922 


1.9480e-07 


-52.365 


6,2, -0.8217* 


0.357 


5.608 


9.2482e-08 


4.551 


6,2, -0.4421* 


1.116 


1.792 


7.8312e-08 


12.779 


6,2,0.0509 


2.102 


0.952 


1.0039e-08 


-83.960 


6,2,0.5463 


3.093 


0.647 


9.0364e-08 


-8.8759 


6,4, -0.464* 


3.072 


1.302 


8.8912e-08 


22.769 


6,4,0.1306 


4.261 


0.939 


2.8022e-08 


-67.972 
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The following parameters were used for the moment of inertia of the ac- 
cretor, the tidal ratio and the mass ratio respectively: I = OAMR^, 
R/d = 0.2 and q = 0.25. It is straightforward to rescale the numbers for 
other values of mass and tidal ratios. The mode quantum numbers marked 
with a * are candidates for tidal synchronization. 



dp5U*dq5U. 



(E6) 



The remaining integrals have been evaluated in iBrinkI l l2005h for 
spheroids of arbitrary ellipticity. In the slow rotation approxima- 
tion, the results are 



a xte ^OpoU OqoU — ZTvp ttt—t-t-, —R 



21 



(fxX^''dpSU*dqSU = 27r/3^ 



1 (Z - m)\ 
{l + -m)\ R 



(E7) 



21 + 1(1- my. w 



X 



2m + (1 - w^){m + l{l + l)w) 



(E8) 



Using l |C9t we arrive at 

MR^n 



4mw H 

w 



- 2 



- 2m + (1 - TO^)(m + l{l + l)w) 



(E9) 



We then obtain u by combining iE9i and l l57t . In table IE II we re- 
port numerical values for A, and iz/Q, for all r-modes with 
m > 0, I < G and I + m even. In evaluating JT and we use 
I — 0AM R^ for the moment of inertia of the accretor, R/d = 0.2 
for the tidal ratio and q = 0.25 for the mass ratio. It is easy to 
rescale the numbers given in table lEl I for other values of /, R/d 
and q. 
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